Unusual features in the nonlinear microwave surface impedance of YBaCuO thin films 
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Striking features have been found in the nonhnear mi- 
crowave (8.0 GHz) surface impedance Zs = Rs+j-Xs of high- 
quality YBaCuO thin films with comparable low power char- 
acteristics (i?res ~35-60 ^D. and Al(15 K) ~ 130-260 nm). 
The surface resistance Rs is found to increase, decrease or re- 
main independent of the microwave field -ff,/ (up to 60 mT) 
at different temperatures and for different samples. However, 
the surface reactance Xn always follows the same functional 
form. Mechanisms which may be responsible for the observed 
variations in Rs and Xs are briefly discussed. 

Measurements of the nonlinear microwave surface 
impedance, Zs, of high-temperature superconductors 
(HTS) is a powerful tool for studying non-equilibrium 
processes in these materials. Nonlinear impedance mea- 
surements allow one to investigate peculiarities of the rf- 
vortex nucleation, and to study the vortex dynamics at 
elevated microwave fields. Such measurements may also 
discriminate between d-wave and s-wave mechanisms of 
pairing symmetry in HTS, and indicate the presence of 
magnetic impurities in the materialsu'D. 

In the present paper, we report observations of non- 
monotonous behavior of Rs and the penetration depth, 
A (or, equivalently, the surface reactance Xs = w/xqA), of 
high-quality epitaxial YBaCuO thin films, in microwave 
fields up to 60 kA/m (^700 Oe) using the coplanar res- 
onator techniquda at 8 GHz. For all samples, depending 
on temperature T, Rs demonstrates completely differ- 
ent behavior, whereas A always preserves the same Hrf- 
dependence, irrespective of sample and T. Measurements 
are presented for very high quality samples over a wide 
temperature range (12-75 K) which at first time reveal 
non-monotonous and uncorrelated behavior in Rs and Xs 
as a function of Hrf ■ Such a behavior does not agree with 
any of thcj-pjijisting models for the nonlinear microwave 
impedanceocl. In the following we discuss several mech- 
anisms relevant to these observations. 

The films are deposited by e-beam co-evaporation onto 
polished (OOl)-orientated MgO single crystal substrates 
10 X 10 mm^. The films are 350 nm thick. The c-axis 
misalignment of the films are typically less than 1%, and 
the dc critical current density Jc at 77 K is around 2 • 
10^ A/cm^. More detailed information on the growth 
technique can be found in Ref.[ |l^. The values of Rs 
and A at 15 K are 60, 35, 50 i^n and 260, 210, 135 nm 
for samples TFl, TF2 and TF3, respectively. 

Changes in Rg and Xs with H^f, ARs = Rs{Hrf) — 
Rs{0) and AXs = Xs{Hrf) - Xs{0), are plotted in fig. |l| 



and fig. ^ for all three samples. For sample TFl for all 
T and in the whole field range ARs ~ H^f, whereas for 
samples TF2 and TF3 the behavior of ARs{Hrf) changes 
dramatically with T. For sample TF2 Rs changes from 
decreasing at 15 K to almost _ffr/-independent behavior 
at 35 K, and to a rapidly increasing function of Hrf be- 
tween 40-75 K. At 15 K ARs diminishes noticeably only 
at Hrf > 10 kA/m showing no features of saturation up 
to the highest available Hrf of ^40 kA/m. At higher T 
(70 K) a transition to a characteristic sublinear field de- 
pendence (^ ^rf^ ri < 1) occurs. A similar behavior is 
also observed for sample TF3 at 15 and 35 K over the 
whole field range, whereas at T > 70 K (see fig. ^) a 
minimum in ARs{Hrf) at low fields appears, after which 
the usual sublinear iJr/-dependence is recovered. As re- 
gards to AXs{Hrf), it is always a sublinear function of 
Hrf at low fields with a characteristic kink and super- 
linear iJr/-dependence (~ H^f,n > 1) at higher fields 
(see fig. |]). This dependence of AXs on Hrf persists for 
all samples and for almost all temperatures, and in gen- 
eral no correlation is observed between AXs{Hrf) and 
ARs{Hrf). The only exception is sample TF3 for which 
AXs{Hrf) qualitatively correlates with ARs{Hrf) at all 
T, and even the minimum at low fields is reproduced in 
both dependences at 75 K (see fig. |l|c and fig. ||c). In 
fig. 1^ and fig. 1^ some of ARs and AXs data are fitted 
to the function ^ H"j which is predicted by Halbritter's 
model of Josephson vortex motion in weak links (WL)i 
(0.5 < n < 2) and the Ginzburg-Landau theory for the 
pair breaking mechanism {n — 2). 

In fig. ^ we plot the temperature dependence of the 
r-parameter (r = ARs/ AXs) for all the samples at three 
microwave power levels. These data are often used-jbo dis- 
tinguish between various nonlinear mechanismsB'H. The 
general trend of r{T) for all the samples is a decrease in 
the absolute value of r with increasing T, gradually sat- 
urating at high T. For samples TF2 and TF3 the most 
pronounced change in r occurs at low T, where it has a 
large negative value of about —1 and rapidly levels off 
with temperature approaching a value of 0.01-0.06 (see 
fig. |b,c). One can see that the initial negative value of 
r is reduced with enhanced power. Unlike other sam- 
ples, TFl in the high field regime {Hrf ~ 7600 A/m) 
displays an increase in the r-parameter at high temper- 
atures (T > 45 K) and reaches a value of ~0.15. In ad- 
dition, the initial low-T r- value for sample TFl depends 
non- monotonously on the microwave field (see fig. ^). 

For explanation of the observed non-monotonous field 
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FIG. 1. Microwave field Hrf dependences of tlie cfiange ir 
tlie surface resistance ARs for tiiree samples TFl, TF2 anc 
TF3 at different temperatures T. The dashed lines are fit 
ting curves using a function ~ -^r/i which is discussed ir 
the text. T- and ri- values are given in the figure. The solic 
line in a) is a fit using the modified model of Gallop et al 
ref[14]. The parameters of the fit are as follows: normai 
resistance R„ = 1.83 f2, zero field critical current density 
Jco = 10^" A/cm^, grain size a — 6.4- 10"^ m, grain penetra- 
tion depth \abo ~ 2.51 • 10~^ m. The insert shows the data 
at 75 K on an expanded scale. 

dependences of ARg and AXg we involve three different 
mechanisms. Each mechanism is capable of describing 
only particular features in 77r/-dependences of Ai?^ aiuji 
AXs- First, the modifiei-weakly-coupled grain modelHl, 
proposed by Gallop et al.t3, assumes that for high-quality 
HTS films a WL between two superconducting grains is 
shunted by another third grain, which serves as an ad- 
ditional path for both dc and rf currents. This model 
presents a highly simplified picture (not least, because 
it makes no distinction between the Meissner and the 
mixed states). Nevertheless it enables one to reproduce 
a reduction of R^, such as we observe, given a certain set 
of the material parameters. However, it is unable to de- 
scribe two important features observed by us; sublinear 
AXs{Hrf) dependence at low fields with a characteristic 
change in curvature at higher Hrf, and the decrease of 
AXs with Hrf (see fig. ||c). We have managed to over- 
come in part the first drawback of the model by intro- 
ducing an effective local flux density Bgff, which inter- 



FIG. 2. The change in the surface reactance AXs as a func- 
tion of Hrf for three samples TFl, TF2 and TF3 at different 
temperatures. The fits are the same as those plotted in fig.l. 
The insert shows the data at 75 K on an expanded scale. 

acts with the rf current. When both the junction and 
grains are in the Meissner or mixed state, the magnetic 
flux is quasi-homogeneously distributed throughout the 
region to which it penetrates, and Beff = B. However, 
in the mixed state of the WL only, the flux will be con- 
centrated inside the junctions due to screening currents 
in the grains, and hence, in the WL B^ff > B. The ratio 
Bef f / B should increase with B and reach a maximum at 
Bf.f f = IJ-qHci {Hci is the lower critical flcld of the grains) 
after which it should decrease rapidly. Adopting a simple 
function for Beff/B{B) which possesses the properties 
specified above (we took the Gaussian function) we have 
managed to get an excellent fit to our AXs{Hrf) data 
(see fig. ||a), but we failed to reproduce the ARg field 
dependence (see fig. |^). Moreover, such a model can 
not reproduce the decrease in AX^ with Hrf observed 
for TF3 sample at high temperatures (see fig. ^). 

Another rrioijiel applicable to our results is the model 
of Eliashberg)13 for superconductivity stimulated by high 
power microwave irradiation. As shown by Eliashberglj, 
in a superconductor with a homogeneous order parame- 
ter distribution, microwave radiation of a certain power 
can induce a new quasiparticle distribution function with 
an increased gap, which in turns leads to an enhance- 
ment in superconducting properties. A similarj-effect is 
predicted by the Aslamazov-Larkin (AL) theorytJ for in- 
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FIG. 3. Temperature dependences of the r-parameter 
(ARs/AXs) for samples TFl, TF2 and TF3 at different H^f 
(specified in the figure). 



homogeneous WL, due to the radiation-induced diffusion 
of quasiparticles out of the junction region which occurs 
for a certain level of microwave power. Since the AL 
theory, contrary to the Eliashberg model, predicts a sup- 
pression of the order parameter at low rf fields, we can 
exclude this mechanism immediately, since we observe a 
reduction of Rg at the lowest fields (see fig. ||3,c). As 
regards to the Eliashberg theory, it predicts a decrease 
of the stimulation effect with lowered temperature, and 
a suppression of superconductivity by a static magnetic 
fieldllj. In fact, for sample TF2 we see that the decrease 
in Rs with Hrf is reduced with increasing temperature 
and completely disappears at high T, whereas for sample 
TF3 the decrease in Rg is observed only at high T (see 
fig. 0c;). Moreover, additional experiments performed by 
us in low dc magnetic fields H^c < 12 mT (to be pub- 
lished elsewhere) , showed that while for sample TF2 iJ^c 
causes an even more pronounced decrease in Rg , for sam- 
ple TF3 the dc field always leads to an enhanced RjiB. 
However, in accordance with Ref . [ [l3[ , stimulation of su- 
perconductivity is expected only in highly uniform nar- 
row and thin superconducting channels with a homoge- 
neous order parameter and microwave field distribution, 
which is hardly the case for our wide and "quasi-bulk" 
samples. 

Finally, the third mechanism which may account for 
our results is the recovery of superconductivity due to 



the field-induced spin alignment of magnetic impurities 
which are likely ta be present in most HTS (particu- 
larly in YBaCuOp. Magnetic impurities are a source 
of Cooper pair breaking due to the spin-fiip scattering 
process. However, at low temperatures the decrease in 
thermal motion leads to the appearance of spin-spin cor- 
relation of the impurity atoms which becomes strong and 
may frustrate the spin-fiip scattering. An external mag- 
netic field also leads to ordering via alignment of the 
magnetic impurity spins and hence, can also lead to a 
reduction of pair breaking. 

Recently Hein et alJ3 observed a correlated reduction 
of Rs and A in both dc and rf fields of the same order of 
magnitude (< 20 mT). They performed an analysis of the 
function AR^/ Rc{AXs/Rc) (where Rc = y/ujfj,o/2an, 
and (T„ is the normal electron conductivity) in terms of 
the two-fluid model (TFM) and found that their data col- 
lapsed onto a single TFM curve. In addition, the conduc- 
tivity ratio y — Gxja^ (where a\ and a-i are quasiparti- 
cle and superfluid conductivities, respectively) was found 
to decrease with increased magnetic field, which was at- 
tributed to the field-induced reduction of pair breaking. 
The majoc, difference between our results and those of 
Hein et alfl is that they did not observe a reduction of Rs 
with Hrj without an accompanying reduction of A; but at 
the same time they observed a reduction of A with Rs be- 
ing almost independent of Hrf ■ In contrast, we observed 
a reduction of Rs for a monotonously increasing A, and 
moreover, only in rare cases did we observe a decrease 
in Rs correlated with a decrease in A (see fig. |l|c and 
fig. ^). In addition, a similar analysis based on the TFM 
was performed by us which showed rather poor scaling 
of our ARs/ Rc{AXs/ Rc) data, as plotted in fig. ^. One 
further distinctive feature of our data compared to those 
of Hein et alB is the significant discrepancy (up to sev- 
eral times) in the Rc values extracted from the ARs and 
AXs data (see table |). Besides, our conductivity ratio 
y = (Ti/(T2, extracted from the fitting to the TFM curve, 
was found to increase with Hrf, rather than decrease, as 
expected for the mechanism of the impurity spins align- 
mentEJ. Nevertheless, our preliminary measurements in 
weak dc magnetic fields (< 12 mT) in field cooled regime 
at constant Hrf showed a decrease of Rs and A with H^c 
for samples TFl and TF2 (to be published elsewhere). 
This suggests that magnetic impurities may play a sig- 
nificant role in our samples and might also affect our 
nonlinear measurements. 

We proceed with an analysis of our high-temperature 

data for increasing|^i?s(i7r/) and AXs{Hrf) in terms 
of the r-parameteii3i3. It is essential to consider not 
only the rjialue, but also the power dependence of the 
impedanceO. A mechanism, such as the response of 
Josephson vortices, for which ARs,AXs ^ H"j (0.5 < 
n < 2) and the r-value is about unity, can be excluded 
immediately, since for our experiments r never exceeds 
0.15. Moreover, the fit of our rf field dependences with 
a function ~ iJ^ (see fig. |^,c and fig. ||b,c) has re- 
vealed uncorrelated values of n for Ai?» and AX, data 
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FIG. 4. Parametric plot of ARs/Rc vs. AXs/Rc for vari- 
ous samples at different temperatures (specified in the figure). 
The solid line is a fit to the two-fluid model discussed in the 
text. The insert shows the low power data on an expanded 
scale. 

(while from the theory^ they should be the same), saying 
nothing about an apparent departure of the fit from 
AXg{Hrf ) data at high fields (fig. Hb,c). The same con- 
clusion is valid for the heating of weak links {rjjE < Ij 
ARs,AXs ^ and the RSJ modefi (rusj < 1, Ai?^ 
increasing and AXs oscillating with i?r/)- A value of the 
r-parameter consistent with our data could follow from 
either uiiiform heating or intrinsic Ginzburg-Landau non- 
linearityO (for both mechanisms r < 10^^), but the rf 
dependence ~ H^^j: is generally not observed for our sam- 
ples (except sample TFI, for which ARg ~ H^j, but AXg 
is not ~ H^j:, see fig. ||a). Moreover, the r- value should 
increase with T for the above two mechanisms, while we 
see almost T-independent behavior at high temperatures 
(fig. ^,c). Thus, our high-temperature data which shows 
an increase in AR^ and AXg with Hrf are apparently not 
explained by any of the known theoretical models. 

In conclusion, in our experiments we appear to observe 
a complicated interplay of several nonlinear mechanisms. 
At low temperatures, the observed reduction in Rg may 
arise due to the effect of the magnetic impurity spins 
alignment by the r/-field, while at higher T stimulation 
of superconductivity by microwave irradiation and vortex 
mechanisms may also come into play. However, universal 
temperature- and sample-independent i?r/-dependence 
of the penetration depth A (or, equivalently, the surface 
reactance Xs) and similar values of r ~0. 01-0. 06 for all 
the samples over a broad temperature range (35 < T < 



TABLE I. Two fluid model fitting parameters of AR^/Rc 
vs. AXs/ Rc dependences for various samples at different tem- 
peratures. Here y{Hrf = 0) and y{H°p) are the conductiv- 
ity ratios {y = ox j 0-2) at low and high microwave powers, 
and Rc{ARs) and Rc{AXs) are the Rc values extracted from 
ARs{Hrf) and AXs(Hrf) data, respectively. 



Sample 


y{H,f = 0) y(//°/) 


Rc{ARs), RciAXs), mO, 


TFI, 15 K 


0.183 


0.173 


0.103 


0.046 


TFI, 35 K 


0.028 


0.123 


0.215 


0.084 


TFI, 60 K 


0.050 


O.Iff 


0.700 


0.352 


TF2, 70 K 


0.008 


0.012 


0.021 


0.004 


TF3, 75 K 


0.010 


0.022 


0.0008 


0.0003 



75 K) do not rule out the possibility that all the observed 
features may arise due to one and the same mechanism, 
the origin of which is not known at the moment. At 
the same time, absence of correlation between ARs{Hrf) 
and AXs{Hrf), for some of the samples particularly TFI, 
implies that microstructure of the samples may interfere 
with the intrinsic behavior in the nonlinear response. 
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